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'OF POCiR ?'’ 

Crop Weather Models of Barley and Spring Wheat Yield 
for Agropbysical Units in 
North Dakota 

Sharon LeDuc 

Previous crop weather models for North Dakota crop reporting districts and 
the entire state were developed for spring \^eat (LeDuc, 1981) and for barley 
(Motha, 1980), This report documents spring wheat yield models for agrophysical 
units in North Dalcota. Agrophyslcal units (APU's) are areas defined to be homo- 
geneous with respect to soils and climate, manifested in crop yield. The defi- 
nition of these areas is contained in a memorandum of understanding from 
Strommen and Dragg (Jan. 10, 198O). Three APU's or groups of counties, were 
defined in North Dakota (Figure 1): west (21), central (19) and east (20). 

Data 

Data were developed for these APU's specifically for use in crop yield 
modeling efforts. Month.ly temperature and precipitation observations for 
cooperative stations (stations within 5 nautical miles averaged and con- 
sidered as a single station) were used to estimate the temperature and precipi- 
tation appropriate for' each county. Counties are the sma3.1 politlca.1 sub- 
divisions in Figure 1. 

The inverse distance method and the Thiessen method (Llnsley, Kohler and 
Paalhus; 1958), were used for averaging the cooperative wee.bher station data. 
Both methods used a weighted average of stations within 50 nautical miles of the 
county's geographic center. This inc3.uded stations in neighboring states but 
not from Canada. The number of stations for a county was zero for a few coun- 
ties in some of the early years but was five for most of the estimates. 
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Figure 1. 

North Dakota Agraphysical Units (APU's) and Counties 

NORTH DAKOTA . * 
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Each station was weighted accordingly in the averaging. The weight (w^) for the 
ith station with the Inverse distance method is w^ « dj[ "I, where d^ is the 
Euclidean distance from the geometric center of the county to that station. The 
Thicjssen method uses weights based on the area closest to that station (Figure 
2). The measure of the area closest to station 3, shown by the lightest shaded 
area in the county (Figure 2 ), is the wei^t assigned to station 3, i.e., W3 = A3. 

This same procedure is used with all three weather variables; maximum 
temperature, minimum temperature and precipitation. Let Vj_ represent one of 
these variables for station i and let station 1,2, . . . , n be all of the sta- 
tions within fifty nautical miles of the center of the county of interest. The 
estimate of this variable for this county is: 

V = k. wiVi/ k. Wl, 
i=l 1=1 


where the weight corresponds to one of the two methods. These data were deve- 
loped under the supervision of Dr. Wayne Decker, Atmospheric Science Department, 
University of Missouri, Columbia. 

The two methods produced similar results. The biggest differences occurred 
for temperature during the winter, due to frontal movements. These differences 
were not of concern since the crops are not planted prior to or during the 
winter. There is no direct physicaJ. explanation that would suggest winter tem- 
perature shoald be included in the nodel. The estimates of weather variables 
from the inverse distance method were used in the nndels because the method is 
much simpler. The counties were assigned to one of the three APU's and the 
value for a single APU was calculated as the simple average of each meteorologi- 
cal variable for each month of each year. These variables formed the basic data 
from which the derived variables were caJ.culated. 
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Selection of Variables 


Predictor variables were selected from the list of potentieQ. predictors 
described earlier. The method of selection involved the use of stepwise 
regression. An attempt was made to examine predictor variables selected for a 
particular APU by the stepwise regression procedure using several selection 
methods: step\'fise, forward and backward. The slgiiflcance level for entry with 

the forward selection was .50, the significance level for remaining in the model 
with backward elimination was .10. Stepwise used the same significance levels. 
The resu3.ts of the three methods were compared and variables included by all 
three methods were considered as potential variables. The rraximum improvement 
method was used and the model selected for which all variables were signifi- 
cant at the .10 level. Several of these models contained a number of related 
variables for the same month. Various combinations were tried in an attempt to 
include information from other months and to reduce the number of variables con- 
taining information from the same month. It is not possible to document here 
all the combinations which were examined nor why many were rejected. The sta- 
tlstica]. significance in the end of season model was a prime consideration. The 
need to include information from different periods in the growing season was 
recognized but not at the expense of Including a statistically Insignificant 
variable. 

The problem of near colllnearity of variables will have to be addressed in 
the eva3.uation of the models. Unless the nature of this near colllnearity 
changes, the precision of the predicted value shou].d not be affected (Mandel, 
1982). Variables significant for adjacent APU's were checked and Included in 
other APU's if statistically significant and physically meaningful. Variables 
Included in the final model were used for the truncated models, i.e., models 


5 


which use only the weatl'ier variables for a stated month and the preceding 
months. These models are used for obtaining estimates within the growin; season 
without making assumptions about future weather conditions. In one case in 
extra predictor variable was Included for May and June truncations but was not 
In the July final model. The variable plus May precipitation minus potential 
evapotransplration was not slgilf leant in the model for the end of the season. 
Variables were considered for Inclusion or elimination as a result of comments 
received from the review of the first draft of this document. 

This was the spring vAieat model for the western APU. This additional variable 
Improves the earlier truncations but do not add slgnlficantl.y to later 
truncations , 

The estifnated coefficients of the squared terms are negative. This is 
intuitively appea!).ing since extrapolation in either direction should at some 
poll- - . associated with a decrease In yields. 

Barley Models 

The statistics defining the models for various truications are shown in 
Table 1. The more Important months for estimating the Impact of weather on 
barley yields appear to be June and July. Evapotransplration appears to be 
a good Indicator of final, yield for June In all APU’s and also in July for the 
Central and Western APU's. Temperature in July is lmport;mt In the Eastern 
May variables did nob seem to be Important and there are no truncated nt> IjIh for 
the month of May. The squared deviation of cumul-ative precipitation (Sept- Apr) 
is important In the Central and Eastern APUhs. The coefficient Is negative 
Indicating that too much or too l.lttle precipitation during tb-' fall, winter and 
spring just prior to planting are two conditions v^ich can be associated with 
decreased yield. 
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Table 1 

Coefficients and for Barley Models (Yield in qu/ha) 
Coefficients for 


Variables 

Area Ihru r 2 TRENDl TREND2 DSQ4 ETP6 BIHAT6 ETHATT DEP7 XDD7 



Trend 

.68 

.507 

*399 





April 

.7b 

7491 

.502 

-.000747 



19 

June 

.83 

.500 

378“ 

-.000496 10.2 



Central 

July 

.94 

.511 

.504 

-.000295 12.9 

-.176 

.0337 


Trend 

.bO 

.359 

35B“ 





April 

.85 

.315 

.620 

-.000354 



20 

June 

.b9 


.624 

-.000311 -.112 



Eastern 

July 

.94 

“351 

7S2b 

-.000227 -.107 


-.916 


Trend 

.51 

.459 

.336 





June 

.75 

.3"C5 

.331 

17.3 



21 

July 

.91 

~743 

358“ 

15.7 

-.109 

1.05b6 

Westem 









TREf>IDl is year minus 1930 up to a maximum of 31. 

TREND2 is year minus I960 down to a mlnirnum of 1. 

DSQ4 is the squared departure of the Sept thru Apr cumulative precipitation from the 
APU average (mm). 

ETP6 is the ratio of actual to potential evapotranspiration for June (unitless). 

ETHAT6, ETHAT7 are the "climatically appropriate" evapotranspirations for June and July, 
respectively (mm). 

DEF7 is precipitation minus the potential evapotranspiration in July (ram). 

XDD7 is the July temperature minus the average (°C). 
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The western APU does not appear to have the problem with excessive 
precipitation. Extremely hot temperatures and moisture stress at heading, 
usual3.y in July, are stron^y related to yield. Similar variables are found in 
the models for the separate APU’s i»iith the central APU showing similarity to 
both the eastern and western APU’s. 

Tile trend variables selected were the same for each model. As mentioned 
earlier, the only trend variables considered were those defined by Motha (198O). 
TRENDl is year minus 1930 until 196I and is 31 after that. TREND2 is one until 
1961 and is year minus i960 after that. 

The coefficients for TREM)1 ranged from .315 to .511 for the various models 
and for TREND2 from .331 to .626. This surrogate for technology indicates the 
contribution from 1948 to 198I for the end of the season models are: 
central APU, + 16.7 quin tals/hec tare; 
eastern APU, + 17.2 qulntalsAiectare; 
western APU, + 12.9 qulntalsAiectare. 

This corresponds to the fact that larger yield increases have been possible 
where nxore moisture is available, i.e., the eastern part of the state. Graphs of 
the fit of the ful.l season models are shown in Figures 3-5. The standard errors 
range from I.60 to 2.49 qAia for the CRD end of season models which were devel- 
oped with data from 1931-1978. Hie standard errors for end of season APU 
models are 1.36, 1.29 and 1.69 qAia for the east, central and west, 
respectively. 

Spring Wheat Models 

The statistics defining the models for various truncations are shown in 
Table 2. The mu3.ticollinearity problem is obvious in the July truncation for 
APU 19 and 20. The coefficient for XDD7 has a large negative coefficient and 
BIHAT7 a large positive coefficient. These effects are offsetting because the 
two variables are highly correlated. 
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Table 2 

Coefficients and r 2 for Spring Wheat MDdels (Yield In qu/ha) 


Area 

Variables 

Thru 

TRENDl TRENDS 

CP4 DSQ4 DEIi5 XDD5 

XDD6 

DEF6 

XDD7 

E3HAT7 defy 


Trend 

.72 

.884 

-.054 







April 

.77 

.915 

.137 

1.000495 





19 

June 

.So“ 

.939 

.134 

7066438 

-.579 




Central July 

.«9 

.903 

.4^4 

1.000386 

-.510 ' 


-678 

134 


Trend 

.74 

T784 








April 

.79 

”7795" 

.453 

-.000324 





20 

June 

rBT" 

■7805” 

7j6I 

-.000292 

-.431 




Eastern July 

' .69 

.763 

77W~ 

7006269 

-.400 


-322 

62 


T?end 

.65 

.732 







21 

April 

.70 

7647 


.296 -.000156 





Western May 


.653 


.0391 -.000549 ".0381 1.10 






June 

.85 

7642 


.0361 -.000407 .0363 0.970 


.0239 




July 

.92 

.684 


.0271 -.000291 766^ 


,0278 

-.46 

.0292 


TREj® 1 is 1 iffitil 1956 , then Is year minus 195“ through 1966 and 12 
TRENDS Is 1 untH 197^, then Is year minus 1972 j 

CP4 Is cumulative precipitation (irm) from previous September thrru^-n April; 

DSQ4 Is the square of the quantity, CP4 minus the average CP4; 

DEP5(DEP6,DEF7) Is precipitation minus potential evapotransplratlon (mm) for May, June, July); 
XDD5, XDD6 (XDD7) Is the temperature minus the average temperature for May (June, July) In ’’C; 
EIEHAT7 Is the climatically appropriate evapotransplratlon for July In ram. 


May truncation la Included .for the western APU only while all three APU's 
have a truncation for April. July> the month of heading, vrould appear to be the 
month when weather has the most Inpact oi yield. The explained variation (R^) 
Increases with the model for the July truncation in comparison to the models for 
the previous truncation, i.e., June. 

The temperature in July is important In the models for all three APU's and 
the coefficient has a negative sign indicating hotter weather In July is asso- 
ciated with decreases In spring wheat yield across the state. The climatically 
appropriate evapotransporation is in the model for the central and eastern APU. 
Ihe western APU Includes the precipitation minus potential evapotranspiration. 
Ihe si@i of the coefficioits for these variables Indicates that more moisture is 
associated v/ith higher yields. 

Cumulative precipitation appears In the models for all APU’s. The cumula- 
tive precipitation from Septent>er prior to planting through April is sqiiared 
after subtracting the average. This variable has a negative coefficient ln.di- 
catlng that too much or too little precipitation during those eight months. 

Which affects soil moisture at planting in April and May, is associated with 
decreased yield. The western APU also has the linear term of the cumulative 
precipitation from September throu^ April. This reflects the fact that this 
APU climatically has less moisture than the other two APU’s, and the impact of 
deficit precipitation is more pronounced, i.e., yields are lower, than with pre- 
cipitation greater than normal. The western APU model also has a variable for 
May which Indicates yield Increases are associated with more precipitation than 
potential evapotranspiration. 

The first trend (TRENDl) for the spring wlneat models was significant for all 
APU’s and is defined as year minus 1954 from 1955 until 1966. It is one 
prior to 1955 and 12 after 1966. An additional trend (referred to as TRENDS) 
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included in the central and eastern APU’s is 1 until 1974, then it is year minus 
1972. Ihe changes in the contribution of trend to the yield as estimated from 
the full season models from 1948 to 1979 are the following: 12.4 quintaie/ 
hectare for central APUi 12.5 quintals/hectare for eastern APUj and 7.5 
quintals/hectare for western APU. Graphs of the fit of the full season models 
are shown in Figures 6-8. 

Comparison with Crop Reporting District (CRD) Models 

Barley 

The trend used in the CRD barley models had two linear components, one from 
1931 until 1961 and the second from 196I through 1978. Five of the nine CRDs 
had an additional quadratic trends from 1961-1978; west central, central, and 
the three southern CRDs. The APU barley ncdels have two linear trends vrtilch by 
df^sign are exactly the same as the linear trends for the CRD models. No quadra^ 
•‘ic trends were included in any of the APU models. None of the APU models have 
May or March tnincations. All except two of the CRDs have models for May 
truncations, and four have a March truncation. Two of the eastern and one of 
the western CRDs have August truncations whereas only the central APU had an 
August truncation. The explained variance (R^) for the full season model is 
very similar, perhaps sli^tly higher (.91 - .94) for the APU models; CRD models 
ranged from .80 to .94. There is quite a bit of difference in the meteorologi- 
cal variables selected as predictors. Ihe July variables did not seem to 
improve the R^ as much with the CRD models as with the APU models and August 
truncations are included for the CRD models. 
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Spring Wheat 

CRD models Included one trend or, in the southeast CRD, tvra trends. The 
trends v^ich were examined for APU's were designed to be the same as in the CRD 
models. The first (IRENDl) Is linear from 1955 until 1966. The second (TREf®2) 
was linear from 1966 until 1973 and constant before and after that period. The 
southeast was the only CRD which included the second trend. An additional 
CIREND3) viiich was considered but not selected by any of the CRD models was 
constant prior to 1973 and linear after that. This trend was selected for the 
eastern and central APU models. All CRDs except the three western and the south 
central had an August truncation. The explained variance R^ in yield for the 
final truncation model, ranged from .85 to .88 for the CRDs, except south central 
which was .77. The explained variances for APU final truncation models were the 
following: .89 central AfU; .89 eastern APU; and .92 western. These were 

quite compatible with the CRD resalts. Although July appeared to be an impor- 
tant month in the APU models, two CRDs (central, east and south central) do not 
even have July truncations and the Ju].y truncation does not improve the 
r 2 appreciably for the southwest CRD. Temperature variaW.es are conuion for 
July, the month of heading, for both types of models. Only the three western 
CRDs and south central CRD have May truncations, correspondingly only the 
western APU model has a May truncation. 
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SUMIIARY 


Ihe mode],s for APUs should provide yield estimates that allow for evaluation 
of the liiprovement over CRD nxjdels. The models are defined in this paper and a 
general comparison is made with the CRD mcdels. A testing procedure will be 
necessary to evaluate the Improvement by using APU Instead of CRD models. 

This study provides the first attenpt to develop agrophysical unit models, 
i.e. models for areas larger than CRDs that are assumed to have similar physical 
characteristics. These candidate models will be evaluated and conpared with the 
CRD models at the state level,. This is a task In the YMD project. Procedures 
for comparing the yield indications at the state level already exist. 
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